Distribution of the area enclosed by a 2D random walk in a disordered medium 
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The asymptotic probability distribution for a Brownian particle wandering in a 2D plane with 
random traps to enclose the algebraic area A by time t is calculated using the instanton technique. 
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It is well known that the properties of random walks 
are dramatically changed by the presence of quenched 
disorder. In particular, if a walker can be irreversibly 
trapped at some randomly distributed sites, then the 
asymptotic probability to return to the starting point (or 
the total probability to survive until time t) is given in d 
dimensions by p(t) ~ exp(— ci d /( d+2 )) [Q, while the mean 
square displacement decreases compared to a pure diffu- 
sion: (r 2 ) ~ t 2 /( d + 2 ) S, On another hand, if a walker 
is advected by a random force, then the total probability 
is conserved, but the diffusion coefficient and the mean 
square displacement acquire logarithmic corrections [3|. 
Further examples can be found, for instance, in Ref. 
Much less, however, is known about the influence of a 
quenched disorder on the topological properties of ran- 
dom walks, which became a subject of theoretical inves- 
tigation very recently ||. By "topological" properties 
we mean such characteristics as the winding number of 
a Brownian particle, or the linking number of a closed 
polymer, or the algebraic area enclosed by the trajectory 
of a random walker. In this Rapid Communication, we 
concentrate on the latter case and calculate the asymp- 
totic distribution of the area swept by a planar random 
walk wandering in the presence of traps. 

The probability distribution for a random walk start- 
ing at some point r' at t = to end at a point r after 
time t obeys the diffusion equation 
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Here D is the diffusion coefficient and U(r) — Uq S(r— 
Hi) is the random "potential" , which represents the trap- 
ping probability per unit time (Uq > 0). The positions 
Ri of traps are distributed uniformly in a plane according 
to the Poisson law with mean density p. The probability 
for a random walk of "length" t to enclose the algebraic 
area A can be obtained by averaging the corresponding 
(5-function constraint over the solutions of Eq. (|l]) in a 
given distribution of traps: 



V{A,t\U) = (8\A- \ I drr\r)e{T) 
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where 9(t) is the angle between the radius-vector r(t) of 
the particle and some fixed direction in the plane. Writ- 
ing the ^-function as an integral over an auxiliary variable 



B, and using the Wiener path-integral formula for the so- 
lution of Eq. (pi), we arrive at the following expression: 
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It is easy to see that the path integral on the right-hand 
side of Eq. (|J) represents the Euclidean Green func- 
tion t/s(r, t; r', 0) of a fictitious quantum particle of mass 
m = (2D) _1 moving in the random potential U(r) and 
in the uniform magnetic field B (we choose the units in 
which h = e = c = 1). Let us now assume that the 
trajectory is closed (r = r') and average the enclosed 
area distribution over all positions of the starting point: 
((...)) r = (I/O) Jd 2 r(...) (Q is the system volume). We 
have: 



T(A,t\U) = dB e lBA Z{t,B\U), 

J —oo 

where Z(t, B) is the partition function at inverse "tem- 
perature" 1/T = t. Averaging now over the positions of 
traps, we finally obtain: 
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where N(E, B) is the average density of states of a quan- 
tum particle described by the Hamiltonian 



H = D{-iV - A(r)) 2 + U(r). 



(4) 



We choose the cylindrical gauge for the vector potential: 
A e = Br/2. 

In an ideal case (i.e. in the absence of traps), the 
eigenvalues of the Hamiltonian (Q) are the Landau lev- 
els E n = u c (n + 1/2), where ui c = 2DB is the cyclotron 
frequency. The density of states is then given by the set 
of equidistant 5-function peaks, and the partition func- 
tion can be easily evaluated, resulting in the following 
expression for the enclosed area distribution: 
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This result was first obtained by Levy quite a while ago 
The form of the dimensionless scaling variable x is 
quite natural, since the only characteristic scale with di- 
mensionality of area in a clean system is the mean square 
displacement (r 2 ) ~ Dt. 

In order to calculate the density of states in the pres- 
ence of disorder, let us first estimate what characteristic 
scales of energy and magnetic field determine the asymp- 
totic behavior of V{A,t). We are interested in calcula- 
tion of the probability distribution at large but fixed t 
and A — > oo. As seen from (^), this limit corresponds to 
E — > 0, B — ► 0, and ui c -C E. Therefore, it looks natural 
to start with the case of B = and make sure that we are 
able to treat the magnetic field as a small perturbation 
in the relevant range of parameters. 

At B = 0, the low-energy behavior of the density 
of states is determined by the rare fluctuations of the 
concentration of impurities creating the large areas free 
of traps which are able to sustain the eigenstates with 
E — ► 0. The asymptotic expression is given by N(E, B — 
0) ~ exp(— const pD/E) at E <C Eq, where E — pU is 
the mean value of the random potential Quantita- 
tively, such exponentially small "tails" are determined 
by the contributions of instantons, i.e. spatially local- 
ized solutions of the saddle-point equations of the effec- 
tive field theory. The typical size of a clean area, or the 
instanton diameter, grows in the time representation as 
Zinst ~ (Dt/p) 1 / 4 [§. If an external magnetic field is im- 
posed on the instanton, then its effect on the energy spec- 
trum can be calculated perturbatively as long as the mag- 
netic length I b = \fl/B considerably exceeds the instan- 
ton dimension k ns t, which is the case if A ^> (Dt/p) 1 / 2 . 
It is this condition that determines the limits of applica- 
bility of our theory. Note also that if we were interested 
in calculation of the intermediate asymptotics of V(A,t) 
at moderate A, then we could use the exact expressions 
for the density of states at E — uj c /2 -C lo c M (since the 
random potential is positive everywhere, the density of 
states vanishes at E < u> c /2). 

Let us now calculate the effect of an external mag- 
netic field on instantons explicitly. The density of states 
can be expressed in terms of the retarded Green function 
of the Schrodinger equation with the Hamiltonian (ji|): 
N{E,B) = -(nfl)- 1 J d D r Im (G R (r, r; E, B))u, where 
G R {r 7 r';E,B) = (r\(E-H + i0)- x |r'). The Green func- 
tion can, in turn, be calculated by standard means of the 
quantum field theory. In order to carry out the averaging 
over the positions of traps, we resort to the supersymme- 
try approach, in which the cancellation of denominators 
is achieved by doubling the degrees of freedom and in- 
troducing the commuting and anticommuting fields on 
equal footing (see, for instance, Ref. ||). Before disorder 
averaging, the retarded Green function can be written as 
the following functional integral: 



G R (r,r';E,B) = -i lim / £> 2 $(r) <p(r)ip(r') 
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where the two-component superficlds 



$(r) 



*(r) = (^*(r),V(r)) 



are composed of a Bose field <p and a Grassmanian field 
if), and X> 2 $ = {l/ir)V(Re<p)V(Imip)V'iijVil>. The prob- 
ability of having N impurities located at the points 
Ri,...,Rat in the area Q in the plane is given by the 
Poisson law: 



Pn(R-i, ...,Rjv) — 
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Using this expression to average (||) over the positions of 
traps, we end up with the following effective action: 



iS[$\ = J d D r |i$ (E - D(-iV - A) 2 + irf) $ 



-P U-e 



)}• 



(7) 



An important property of the action (Q) is that it is in- 
variant under the global supersymmetry transformations, 
mixing the boson and fermion sectors: 



(8) 



where T is a 2 x 2 unitary supermatrix [Bj. Due to this 
symmetry, the problem of finding the saddle points of (Q) 
can be considerably simplified. Indeed, since the saddle 
point manifold is invariant under the transformations (|^), 
we are able to seek the instanton solution in the following 
form: 
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where ^inst(r) = <^mst(0 is a cylindrically symmetric spa- 
tially localized function. If to write the boson fields as 
(p = tpi + i(p2, <p* = <pi — iip2, where <p\p. are real on the 
initial functional integration contour, then the imaginary 
part of the Green function is determined by a non-trivial 
saddle point of the action in the complex plane of tp± 2 
0,0, and 
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with the exponential accuracy. The fermion sector can be 
neglected as long as we are not interested in calculation 
of the pre-exponential factor in (|l(]) . 

In order to make Eq. (^) in the boson sector real, we 
rotate the integration contour: tp\^ — > 



t/4 



As 



result of this, the exponent on the right-hand side of Eq. 
(0) changes: iS[Q>] — > — where the action S is 
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+p(l-e-^)} (11) 
(i = 1,2). Introducing the dimcnsionless variables: 
r = £x, <Pinet(r) = y/Uof(x), 



where £ = yjD/E, we obtain, from Eq. (|TT[) , a non-linear 
differential equation for the saddle-point solution /: 

x dx \ ax J 



where a = yjE a jE > 1 and = uj c /AE < 1. Similar 
equations without magnetic field were obtained in Refs. 
J[^ , [L3| , using different techniques. 

Due to complexity of the equation (|l2|), we are able 
to find only an approximate solution. To this end, we 
replace the "potential" V(f) = a 



— n/ 2 p f 



in Eq. (Dj) by a 



piecewise constant potential: 



Ves(f) 



at / < 1, 
at / > 1. 
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Then, Eq. (|12|) becomes effectively linear and reduces to 
a couple of the Schrodinger equations, whose solutions 
f 1,2(2) satisfy the following conditions: 

fi(xi) = f 2 (xi) = 1, fi( Xl - 0) = f' 2 (x x + 0), (14) 

where x\(a,0) is the position of the discontinuity in 
the effective potential, which is to be determined self- 
consistently. Going back to the dimensional variables, it 
is easy to see from ([ll]) and (|l2|) that the instanton ac- 
tion is proportional to the area of the effective potential 
well: 



irp^xj 



(15) 



Let us start with the case of B = 0, i.e. (3 = 0. The 
solution of the linearized saddle-point equations can be 
written in form 

f m = { AW = C i J oW, at0<x<xi, , > 
K ' \ f%(x) = C 2 K (ax), at x x < x, K > 

where Jo(x) and Kq(x) are the Bessel functions of real 
and imaginary arguments, respectively. Substituting this 
solution in the matching conditions (|l4|), we obtain the 
equation 



Jx(x) K 1 (ax) 



a- 



J (x) K (ax) 
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In the limit a — > 00, assuming that xi *~ 1 and using 
the asymptotic expansions of the Bessel functions juj, 
we obtain Jq(xi) = 0, i.e. x\ = a, where a ~ 2.405 is the 
first zero of the function Jo (a;). After substitution in jl5| ) 



and (10), the Lifshitz result N(E) ~ exp(— const pD/E) 
is recovered. 

In principle, one can find the exact solutions of the 
Schrodinger equations inside and outside the potential 
well at B ^ (they are expressed in terms of the conflu- 
ent hypergeometric functions), match them at the point 
x = X\ and finally end up with a transcendental equa- 
tion for x\ (a, 0), which can be solved at — > 0. However, 
we prefer not to follow this procedure here, because the 
same results can be obtained using much more physically 
apparent reasoning in the spirit of the Lifshitz 's original 
derivation 0. First, we note that, in the limit a — > 00, 
the instanton solution satisfies the Schrodinger equation 
for a particle confined in the potential well with infinitely 
high walls in a magnetic field. The ground state energy 
eo is equal to unity (in the units of E) . In the absence of 
magnetic field this condition fixes the radius of the well 
at xi(0) = a. At 7^ 0, the lowest order perturbative 
correction to the ground state energy is 



2 I dxx 3 f(x) 
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where f(x) ~ Jq( x ) is the unperturbed ground state 
wave function ([lj]). Calculating the integrals with the 
Bessel functions, we obtain: 6cq = c/3 2 , where c ~ 1.261. 
To keep the ground state energy fixed, this correction 
should be compensated by the corresponding increase in 
the radius of the potential well: x\(0) = x\{Q)(l + 5eo). 
Substituting this in (|l5|), we finally obtain: 



Sinst{E 7 B) 



irpDa 2 
E 



cD 2 
4E 2 



B' 



(19) 



This expression is valid atE'^O, B — > 0, w c /E —> 0. 

The asymptotic probability distribution of the enclosed 
area can now be obtained from (||) , ([l(]) and ( |l^ ) by cal- 
culating the integrals by the steepest descent method: 

/CO POO 
dB / dE e iBA e -Et e -S iaet (E,B) 
-00 JO 

-exp^v/^^expj-^^}. (20) 

The first exponential on the right-hand side is noth- 
ing but the asymptotic "tail" of the total probability 
p(t) — fdAV(A,t) for a random walker not to be 
trapped M. The second exponential thus represents the 
conditional probability to enclose the area A, provided a 
walker has survived until time t. 

We see that the asymptotic behavior of V(A, t) is dras- 
tically changed by the presence of disorder, compared to 
the Levy's result (ph. The distribution becomes Gaussian 
V(A,t) ~ exp(— x ), with the scaling variable 



(21) 
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so that the standard deviation now grows slower than in 
the clean case: (A 2 ) 1 / 2 ~ t 3 ^ 4 (the mean value (A) is, of 
course, zero). Such a different form of the scaling variable 
can be related to the presence of an extra length scale 
r p ~ y^Jp hi the system, which depends on the concen- 
tration of traps (nothing depends on the absolute value 
Uq of the random potential). Note also that a similar 
Gaussian distribution was obtained in Ref. jl5| for a 2D 
random walk in a box of a finite size L. In this case, the 
fictitious magnetic field was also treated perturbatively, 
giving rise to (A 2 ) 1 / 2 ~ L(Dt) 1 / 2 . Qualitatively, such a 
difference between the two systems is due to the fact that 
in our case the size of the effective potential well is not 
constant, but grows with time as Lit) = k ns t(t) ~ t 1 / 4 . 

In conclusion, we studied the asymptotic probability 
distribution of the algebraic area enclosed by a planar 
random walk in the presence of immobile random traps. 
It is shown that this probability is directly related to 
"the Lifshitz tail" in the density of states of a quantum 
particle in a Poisson disorder and uniform magnetic field. 
In contrast to the case of an ideal random walk, the en- 
closed area distribution turns out to be Gaussian with 
the standard deviation growing as 
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